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In this paper we introduce a design for an optical topological cluster state computer constructed 
exclusively from a single quantum component. Unlike previous efforts we eliminate the need for 
on demand, high fidelity photon sources and detectors and replace them with the same device 
utilised to create photon/photon entanglement. This introduces highly probabilistic elements into 
the optical architecture while maintaining complete specificity of the structure and operation for 
a large scale computer. Photons in this system are continually recycled back into the preparation 
network, allowing for a arbitrarily deep 3D cluster to be prepared using a comparatively small 
number of photonic qubits and consequently the elimination of high frequency, deterministic photon 
sources. 



I. INTRODUCTION 

Designing and constructing a viable large scale quan- 
tum information processing system has been the focus 
of extensive research ever since the first basic architec- 
tures were proposed in the 1990's [TMlOj. Experimental 
progress towards this goal has been pronounced [TTHT9] 
and researchers are becoming slowly optimistic about the 
future potential for building a large scale computer. 

In recent years, the extraordinary advances in exper- 
imental systems and theoretical techniques for quan- 
tum information processing has allowed for serious ques- 
tions in architectural design and construction to be dis- 
cussed [2"0H2"T] . This new area of research is generally be- 
ing referred to as Quantum Engineering (QE). Broadly, 
the primary goal of QE is to adapt and combine the best 
experimental technologies and theoretical techniques to 
construct an experimentally viable large-scale quantum 
computer. 

One proposed architecture was introduced in 2009 20J . 
This optics based computer is based on the topological 
cluster state model of computation 28 'M)[ and a photon- 
atom-photon coupling device called the photonic mod- 
ule [3T]. This architecture illustrated the structure and 
operation of a fault-tolerant, fully error corrected quan- 
tum architecture. However, the architecture was based 
on components which were all theoretically determinis- 
tic. Deterministic photon-photon coupling utilized the 
non-linearity afforded by an atom/cavity system present 
within each photonic module and photonic sources and 
detectors were simply assumed to be deterministic and 
of high fidelity. 

It has been well known since the seminal paper of Knill, 
Laflamme and Milburn [BJ that an all optical quantum 
computer based on linear elements only allows for cou- 
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pling between qubits in a probabilistic fashion [321 133] . 
Since this result there has been extensive work inves- 
tigating how probabilistic techniques could be utilized 
to construct a viable architecture [TO l [34ti44] , However, 
while this work demonstrated that in principal proba- 
bilistic components could be used to slowly grow large 
entangled states suitable for quantum computation two 
problems remained. The first is that two dimensional 
cluster states |45j do not include any protocols for quan- 
tum error correction. This problem was addressed by ap- 
plying error correction protocols on top of the underlying 
cluster model and resulted in extremely high demands on 
quantum resources [40] . The second significant problem 
is that these results did not show how such a massive 
optical system was to be constructed, arranged and op- 
erated when large scale algorithms and error correction 
require billions of photons and millions of optical compo- 
nents. 

The experimental development of optical systems |19[ 
I46H52] (and more generally, probabilistic quantum com- 
ponents [S31 [SI]) has in many ways been far more suc- 
cessful than deterministic technologies [T31 [T3J [SSHSH]- 
Therefore, an important problem is can these more ad- 
vanced probabilistic technologies still be integrated and 
used in viable large scale quantum information architec- 
tures. The majority of research attempts to do this from 
the bottom up, taking well established non-deterministic 
protocols and incorporating them into appropriate com- 
putational and error correction models [BUI EI]- While 
this addresses some issues related to more effective use 
of error correction techniques, these results still do not 
address the fundamental architectural structure of a large 
scale system. 

We will attempt to approach this from the top down. 
The optical architecture proposed in Ref. 20 is designed 
such that its physical structure and operation can be very 
well defined up to billions of qubits. Since the archi- 
tectural structure of this system is so well defined, this 
will be our starting point. The system is constructed 
from three key deterministic components. We will show 
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how we can remove two of them. We will take the opti- 
cal architecture consisting of deterministic single photon 
sources, deterministic coupling via the photonic module 
and deterministic single photon detection and replace the 
sources and detectors with the photonic module. This 
will lead to a network running with highly probabilistic 
single photon sources and entirely constructed from one 
quantum component, namely the photonic module. 

To combat the issue of probabilistic sources we intro- 
duce a perpetual design. As the photonic module can act 
as a non-demolition photon detector, photons are simply 
recycled. In this way, probabilistic sources are responsi- 
ble for two tasks. (1) Providing the photons to initialize 
the network and (2) replacing photons which are period- 
ically lost during computation. This paper will demon- 
strate how a highly probabilistic source can be integrated 
into a large scale architecture without sacrificing perfor- 
mance or the overall design and operation of the system. 
This represents the first step at integrating probabilis- 
tic technologies with deterministic technologies and will 
hopefully lead to future architectural designs incorporat- 
ing an increasing amount of probabilistic components. 

We begin in section [TT] with a brief discussion of the 
general principals governing the topological cluster state 



cal qubit, initially prepared in the |+) = (|0) + |l))/\/2 
state, and each edge represents a controlled-cr z entan- 
gling gate between qubits. This unit cell of the cluster 
extends in all three dimensions, dictating the size and 
error correcting strength of the computer. Computation 



model. In section 
troduced in Ref. 



III| we review the optical computer in- 
20j . Section |IV| discusses how we can 
reduce the required quantum technologies by replacing 
single photon sources and detectors with the same device 
utilized in the preparation of the cluster. Section|V]intro- 
duces the idea of a perpetual architecture design, where 
a comparatively small number of photons arc recycled 
again and again to perform large cluster computations. 
Each element of the architecture is detailed, including 
how fault-tolerance is maintained and how photon loss 
is reliably detected. The network consisting of photonic 
sources is discussed in Section |VI| where highly proba- 
bilistic sources are used to "boot-up" the computer and 
replace heralded loss events during computation. Finally 
in sections. |VII| and | VIII| we present several network sim- 
ulations illustrating that a perpetual computer utilizing 
photon recycling and highly probabilistic sources can op- 
erate effectively. 



II. TOPOLOGICAL CLUSTER STATE 
COMPUTING 

Topological cluster states were introduced by 
Raussendorf, Harrington and Goyal in 2007 [28l l29] . 
This model incorporates the ideas from topological quan- 
tum computing [55] and cluster state computation [15] , 
This model for quantum information processing is very 
attractive as it incorporates fault-tolerant quantum 
error correction by construction and exhibits a high 
fault-tolerant threshold. 

Computation proceeds via the initial construction of 
a highly entangled multi-qubit state. Fig. [T] illustrates 
a unit cell of the cluster. Each node represents a physi- 




FIG. 1: Unit cell of the 3D cluster required for topological 
cluster computing. Each node represents a physical qubit ini- 
tialized in the |+) = (|0) + |l})/\/2 state and each edge is 
a controlled-(j z operation between two qubits. The 2D cross 
section of the cluster dictates the size (and error correcting 
power) of the computer while the third dimension dictates the 
total number of computational time steps available. Compu- 
tation proceeds by measuring the cluster along one of the 
three dimensions. 

under this model is achieved via the consumption of the 
cluster along one of the three spatial dimensions [25J 130] 
(referred to as simulated time). Logical qubits are de- 
fined via the creation of "holes" or "defects" within the 
global lattice and multi-qubit operations are achieved via 
braiding operations (movement of defects around one an- 
other) as the cluster is consumed along the direction of 
simulated time. Fig. [2] illustrates a braided CNOT oper- 
ation. Qubits within the cluster are selectively measured 
in the <j z basis to create and manipulate defects. By mea- 
suring the correct set of physical qubits, defects can be 
moved as the cluster is consumed. Physical qubits not 
associated with defects are measured in the a x basis and 
are utilized to perform fault-tolerant error correction on 
the system. 

The specific details for computation under this model 
are not important for this discussion and we encourage 
the reader to refer to Refs. [551130] for further details. For 
this analysis, we will be examining the network required 
to successfully create the entangled cluster for computa- 
tion. 



III. OPTICAL TOPOLOGICAL COMPUTER 

Shown in Fig. [3] is the basic structure of the prepara- 
tion network for the computer. The preparation network 
consists of two sets of fabricated "wafers" containing an 
interlaced network of photonic modules [31] oriented at 
90° to each other [Fig. [3k)]. Each set of parallel wafers are 



FIG. 2: Diagram of a logical CNOT operation in the topo- 
logical cluster model. Each point in the diagram represents 
a physical qubit and the cluster is consumed from the front 
of the image to the rear. This image shows a CNOT which 
is approximately 50% complete. Logical qubits are defined 
via pairs of defects. Four sets of defects are required to per- 
form the braided CNOT. The three blue defects represent the 
control input, control output and target respectively. These 
three defects are known as primal defects |29II30| . The purple 
defect is known as a dual defect and is used to enact braided 
logic operations (logic operations can only be performed be- 
tween defects of differing type). The total size and separation 
of each defects dictates the error correcting power of the topo- 
logical code. 



not interconnected, the only connections are at the junc- 
tion between horizontally and vertically oriented wafters. 
In total there are four separate stages in the preparation 
network, two on each horizontal wafer and two on each 
vertical wafer. The computer operates via the injection 
of single photons into the left hand side of this network. 
Each photon then interacts with a total of four individual 
photonic modules which acts to deterministically entan- 
gle photons into the required cluster. After photons are 
entangled, they are then measured in appropriate bases 
to perform fault-tolerant, error corrected quantum com- 
putation. 

Shown in Fig. [3]d) is a schematic of one of the wafers. 
This wafer illustrates stages three and four of the cluster 
preparation network and an additional array of photonic 
modules that are utilized to perform measurement of the 
cluster (this will be discussed in Sec. VB|. The total 



width of each of these wafters remains constant, while 
the length of each of these wafters scales linearly with 
the size of the computer. 

The injection of single photons into the computer re- 
quires a specific temporal arrangement. Illustrated in 
Fig. [4] is the temporal arrangement of photons for stages 
one and two of the preparation network for one of the 
vertically oriented wafers shown in Fig. [3^l) . Within the 
network, half of the optical lines contain photons tempo- 
rally separated by 2T and the other half have photons 
temporally separated by 4T, where T is the operational 
time of a single photonic module. This arrangement al- 




FIG. 3: Structure of the optical topological computer intro- 
duced in Ref. |20| . The optical computer is constructed via a 
stacked array of "wafers" . Each wafer consists of a network 
of photonic modules 31 , each of which is designed to deter- 
ministically couple single photons. Subfigure a) illustrates the 
structure of the actual computer, two sets of wafters, oriented 
at 90° to each other are connected to prepare the required en- 
tangled cluster. Photons travel through the network from left 
to right through four stages of cluster preparation before they 
are measured. Subfigure b) illustrates stages three and four 
of the cluster preparation network and a second array of pho- 
tonic modules used for measurement [Sec. V B| . The width of 
each wafer is independent of the total size of the computer 
while the length increases linearly with the size and/or error 
correcting power of the computer. 



lows the photon stream to essentially "flow" through the 
interlaced network of modules, maintaining temporal sep- 
aration and ensuring that only a single photon is present 
within any photonic module at any given time. Utiliz- 
ing this temporal arrangement, it was shown in Ref. |20j 
that this network could deterministically prepare an ar- 
bitrarily large 3D topological cluster without employing 
sophisticated photonic routing or storage. Additionally 
the specific operation of every photonic module in the 
network is completely specified and independent of the 
size of the computer. 
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FIG. 4: Temporal arrangement for stages one and two of the cluster preparation network. Half of the optical lines have photons 
separated by 4T while the other half of the optical lines have photons separated by 2T, where T is the operational time for 
a photonic module [3T|. This temporal arrangement allows for each photon to flow through the network of modules. After 
exiting the network, each photon is immediately measured to perform computation. 



IV. REDUCING REQUIRED TECHNOLOGIES 

The optical computer requires several different, high 
fidelity, quantum components. (1) The Photonic Mod- 
ule [31] . (2) High frequency, on demand, single photon 
sources. (3) High fidelity single photon detectors. (4) 
High fidelity single photon switching. It could be argued 
that the construction of large numbers of photonic mod- 
ules represents the most difficult element of the above list. 
If the construction of such a device could be achieved, it 
could be safely assumed that the other required quan- 
tum components will also exist with sufficient reliability 
to construct a large scale computer. However, what if 
this is not the case? Can we further reduce the required 
technology and still design a viable architecture? 

Aside from the photonic module, the next two compo- 
nents listed above that will require extensive development 
are high frequency, on demand, single photon sources and 
high fidelity single photon detectors. These components 
currently exist at several levels of efficiency and relia- 
bility [TQl . However, these technologies still require 
significant work before they are adaptable to a large scale 
quantum computer. These are the two components that 
we will eliminate. We will replace high fidelity single pho- 
ton detectors with the photonic module itself and we will 
replace the high frequency, on demand, sources with non- 
deterministic sources. While we could consider multiple 
non-deterministic sources, we will assume that photonic 



modules are used, distilling weak coherent light into sin- 
gle photons with a very low probability of success. In 
this way we will redesign the computer architecture to 
consist of only two required components, a high fidelity 
photonic module and reliable single photon switching. 

It should be stressed that while the goal of this re- 
design may seem theoretically trivial, from an architec- 
tural standpoint it is significantly more complicated. In 
principal, any number of techniques could be utilized to 
remove the requirements of single photon sources and de- 
tectors. However, we need to accomplish this under some 
tight constraints. 

1. Original design of the architecture required one 
high frequency (MHz or above) source per opti- 
cal line or an ultrahigh frequency source (THz or 
above) serving multiple optical lines. When replac- 
ing deterministic sources with non-deterministic 
sources, this ratio of the number of sources per op- 
tical line must remain effectively constant. 

2. Non-deterministic sources such as the photonic 
module prepare single photons at random (but her- 
alded) times. As these photons may be required 
anywhere in the network, the design must allow 
us for this without introducing complicated photon 
routing or significant photon storage. 

3. The probability of photon loss within the computer 
is non-zero. Therefore the design must be able to 



5 



effectively replace lost photons, which will occur in 
random but heralded locations, with new photons 
which will also be prepared at random but heralded 
locations. 



A. Utilizing the photonic module as a detector 

The ability to utilize the photonic module to perform 
individual photon measurement is the key to replacing on 
demand sources without reintroducing additional tech- 
nologies. As described in Ref. [M] the sole action of 
the photonic module is to project an arbitrary TV-photon 
state, jy into a ±1 eigenstate of the operator X® N . 

Within the preparation network of the computer, each 
module interacts with five photons before being measured 
and decoupled. This, combined with suitable local opera- 
tions, projects each five photon group into ±1 eigenstates 
of the operator ZZXZZ. These are the relevant eigen- 
operators describing the topological cluster state utilized 
for computation [5U1 [Ml US] ■ If, however, only a single 
photon is allowed to interact with a photonic module be- 
tween its initialization and measurement, an arbitrary 
state \tp) is mapped to, 

Mm+) a = l -m+xm\gi)a + \m - *iv»i<?2> a 

(i) 

where M is the interaction between the atomic system 
in the module and the photon, |<?i) a , |<?2)a are the two 
states of the atomic qubit and |+) a = (|<7i) a + |<?2)a)/v2- 
Measuring the atomic system will project the photon into 
either the |+) = (|0) + |1»/V2 or |-> = (|0) - |1»/V2 
state dependent on measuring the module in the \g\) a or 
\g2)a state respectively. Therefore, the module can be 
utilized to perform an X basis measurement on any pho- 
ton it is allowed to interact with between initialization 
and measurement. Combining this with appropriate lo- 
cal rotations before measurement (via optical waveplates 
or waveguide techniques |47j ) allows for the measurement 
of a single photon in any desired basis. 

The major advantage to using the module is that it 
is intrinsically a non-demolition measurement. The pho- 
tonic state is measured via readout of the atomic sys- 
tem. Therefore the photon is not physically destroyed 
during measurement and will simply exit the module in 
exactly the same way as for the preparation network. 
Consequently, we can recycle it. Measured photons can 
therefore be rerouted back into the input of the cluster 
preparation network. 



these photons still need to be initially prepared and loss 
events within the network need to be compensated with 
some type of source device. 

In Ref. [30] it was assumed that the cluster network 
was fed via high frequency, on demand, single photon 
sources. Assuming one such source per optical line in 
the network, the operational frequency of each source 
between 500KHz and 100MHz (assuming that the pho- 
tonic modules have an operational time between 10ns and 
1/xs [3TJ[64]). Replacing a high frequency source with a 
probabilistic source and photon recycling can, in many 
circumstances, be desirable. In the case of our architec- 
ture, the photonic module can act as this probabilistic 
source. This allows for a large scale quantum architec- 
ture to be constructed using essentially only one quantum 
component. 

The photonic module acts as a parity check device, al- 
lowing us to effectively measure multi-qubit observables. 
It is well known that such a quantum component can be 
utilized for distilling weak coherent light into single pho- 
ton states |65|I66|. A weak coherent pulse, in the number 
basis, can be approximated as, 



|0)+a|l) + — 12) + — 13) (2) 
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for a <C 1. This state can routinely be prepared in the 
laboratory. 
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FIG. 5: Internal structure of the atomic system utilized in 
the photonic module. The two ground states, \gi) a and \g%) a 
are the two levels measured when each module is read-out. 
The transition between \gi) a and a third state \e) a is detuned 
from the cavity mode by an amount A. Utilizing a basic 
Jaynes-Cummings interaction in the dispersive limit, a phase 
shift will accumulate on the state \gi) a if the cavity mode is 
occupied with a single photon. 



B. Using the photonic module as a probabilistic 
source 

By using photon recycling we can, in principal, oper- 
ate a large topological cluster state computer with a com- 
paratively small number of individual photons. However, 



The internal mechanism of the photonic module is 
based on working in the dispersive limit of the Jaynes- 
Cummings model. Illustrated in Fig. [5] is the structure 
for the internal atomic system for each photonic module. 
In the dispersive limit, the effective Hamiltonian of the 
system is, H — (3a'ao~ z . The operator(s) a (a') are the 
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creation (annihilation) operators for the cavity mode, <j z 
is the usual Pauli operator over the qubit spanned by 
the two ground states \gi) a and \g2ja an d the effective 
coupling constant is ft = —g 2 /A. If the cavity mode is 
detuned by A from the transition |<?i) a ■<-> \e) a and the 
atomic system is initialized in the state (|<7i) a +|S2)a) / v2j 
then a single photon introduced into the cavity mode will 
result in the following evolution, 



1 

71 



(|fll). + |«&)a)|l> 



1 

71 



2i/3t 



gi)a + \g2)a) |1), 



.(3) 

for an interaction time t. If we tune the interaction time 
to t = 7r/(2/3) before removing the photon from the cavity 
and out-coupling it to a waveguide, we will induce a tt 
phase shift on the state \gi) a - 

Instead of a single photon we now interact the atomic 
system with a weak coherent pulse. The system evolution 
is, 



-y= (|fll)a + l52>a) l«) ^^Jf 



ae 2 



V2 



(|fll)a + \g 2 )a) |0) 

l52>a)|l> 



2 ~' 

are 2 



(|<7l)a + \g 2 )a) |2) 



" e 7^(-|.9l>a+l.92>a)|3>. 



Defining the states, |±) Q 
rewriting Eq. [4j we have, 



2v/3 

(4) 

= {\g 2 )a ± |5i)a)/V2, and 



e = |+)«(|0) + ^|2)) + oe=V-|-) (|l) + -^|3)) (5) 

The atomic system is measured in the |±)„ basis. The 
probability of each outcome is |68) , 



Therefore, if the atomic system within the module is 
measured in the |— ) a state, the projected optical state 
approximates a single photon. The strength of the weak 
coherent state is the determining factor in how close the 
projected eigenstate is to a single photon, with a fidelity 
given by, 



F=|(lhA)| 2 = 
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(8) 



The relationship is inverted when considering the prob- 
ability of measuring the module in the |— )„ state. As 
a — > 0, the probability of measuring the module in the 
|— )o state and projecting the coherent state into an ap- 
proximate single photon state approaches zero. Hence, 
there is a tradeoff between the probability that a module 
will successfully distill a single photon and how well the 
distilled state approximates a single photon. 

Utilizing the photonic module to distill weak coher- 
ent states will therefore result in sources with very low 
probabilities of success (as we require distilled states to 
approximate single photons to a high degree). However, 
even with low success probabilities, modules combined 
with photon recycling can be combined successfully. 



V. PERPETUAL NETWORK 

The previous section illustrated how the photonic mod- 
ule can be utilized to effectively replace single photon de- 
tection and probabilistically distill weak coherent states 
into single photons. We can now discuss the general 
structure of a perpetual quantum computer. Fig. [6] il- 
lustrates the overall structure of the design. Probabilis- 
tic sources are used to slowly "boot-up" the computer. 
Each injected photon then proceeds through the prepara- 
tion network of photonic modules, through a network of 
detection modules and then rerouted back to the source. 
We examine each component of the network separately, 
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and the resultant states after each measurement result is, 
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A. Preparation Network 

The preparation network has already been detailed in 

(6) Ref. gDJ. The temporal arrangement of photons is such 
that half of the optical lines consist of photons separated 
by 2T and the other optical lines contain photons sepa- 
rated by AT, where T is the operational time of the mod- 
ule. Each photon interacts with four separate modules 
and suffers a delay of T for each interaction. 

Within the network there are temporal "windows" 
which allows for each module to be measured and reini- 
tialized. Both the measurement window and the reinitial- 
ization window arc assumed to also take time T. After 
the preparation network, there is an additional delay of 
AT before each photon enters the measurement network. 
This delay allows the final parity checks in cluster prepa- 

(7) ration to be completed before photons are measured. 
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FIG. 6: General structure of the preparation network. Ded- 
icated single photon sources and detectors are replaced with 
photonic modules. As each module acts as a non-demolition 
detector, photons can be rerouted back into the start of the 
preparation network. Each deterministic source can then be 
replaced with low probability sources which provides the ini- 
tial photons to saturate the preparation network and replace 
photons lost during computation. 

B. Detection Network 

As shown in Sec. |IV A| the modules can be utilized to 
perform non-demolition detection on each photon. How- 
ever, there is issues related to fault-tolerance that needs 
to be addressed when designing the detection system. 

It is well known that photon loss is a major error chan- 
nel for optical computers. The topological cluster codes 
that are used in this computer are quite efficient at cor- 
recting this type of error. Several recent results have 
examined the robustness of the topological cluster model 
when subjected to loss [57], suggesting that this partic- 
ular error channel is preferable over standard quantum 
errors. 

In the original design of the architecture, photon de- 
tection was achieved using dedicated single photon detec- 
tors which destructively measure photons and hence can 
discriminate between the presence or absence of the phys- 
ical photon. This combined with the fact that new pho- 
tons are continuously injected from deterministic sources 
meant that loss was an easily correctable error for the 
original design. However, moving to a perpetual archi- 
tecture presents two problems. 

1. The internal function of the photonic module is 
such that it cannot discriminate between a photon 
in the |+) Q state and the vacuum. 

2. As photons are being recycled, undetected loss er- 
rors will be temporally correlated during computa- 
tion (since the same physical photon is continuously 
re-entanged into the cluster at later times). 

Therefore, we need to first figure out a method to 
uniquely detect loss events using the modules and ensure 



that this technique can be made fault-tolerant (i.e. not 
cause errors in detection to spread to large groups of er- 
rors at later times). The detection network is illustrated 
in Fig. 
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FIG. 7: Detection network of photonic modules. Illustrated 
is a small cross section of the network. The upper and lower 
optical lines run at a temporal separation of 4T while the 
central line operates at a temporal separation of 2T. Loss can 
be uniquely identified by measuring photons twice separated 
by a single photon phase rotation. As each module requires a 
temporal window of T for measurement and reinitialization, 
the central optical line has two sets of measurement modules. 
While one module is reinitialized, the next incident photon is 
routed to a second module that was initialized in the previous 
time step. 

Illustrated is a small cross section of the network. The 
upper and lower optical lines run at a repetition rate of 
4T while the central optical line runs at a repetition of 
2T. As we assume each photonic module requires a win- 
dow of T for both measurement and reinitialization, the 
detection network for the central optical line has twice 
as many modules. While one of the modules is being 
reinitialized for the next measurement, the next incident 
photon arrives at the same time. Therefore, it is switched 
to a second module which has already been initialized. 

Fault-tolerance is achieved via measuring photons 
twice. We can write down the module transformations 
for an incident photon in the |+) state, the |— ) state and 
the vacuum. 

M\+)\+) a = \+)\+) a 

M|-)|+) = |-)|-) (9) 

M|vac)|+) a = |vac)|+) a . 

Hence, measuring the atomic system in the |+) state 
indicates the presence of either a single photon |+) state 
or the | vac) state. 



In order to discriminate between the states |+) and 
| vac), we measure the photon a second time. Before 
this second measurement a phase rotation is applied to 
the photon, taking |±) 4-> |=f). Therefore, if the state 
entering the measurement network is the vacuum, both 
modules will be measured in the |+) state. Assuming 
that measurements of the photonic module are error free, 
we are therefore able to uniquely discriminate the actual 
photonic states, |±) and the loss channel |vac). 

Although the above scheme allows us to uniquely iden- 
tify loss in the network, we also need to check that it is 
still effective when each of the two photonic modules are 
subjected to measurement errors and when loss can oc- 
cur between the two measurements. We can summarize 
the possible measurement outcomes for the two modules 
under these error channels. 



Scenario 


Results 


|vac) 

|+) & error on M2 
\ — } & error on Mi 
+} & loss after Mi 
|— ) & error on Mi & loss after Mi 


Mi +) Mo 

+)m 2 
+)m 2 
+)mi +)m 2 
|+)Mi +)m 2 


1+) 

-} & error on Mi & M 2 
|vac) & error on M2 
|+) & error on Mi & loss after Mi 
-) & errors on Mi & M 2 & loss after Mi 


+)mi|-)m 2 

+)mj — )m 2 
+)mi -)m 2 

|+)Mi|-)m 2 
+)Mi|-)m 2 


R 

j+) & error on Mi & M 2 
|vac) & error on Mi 
|+) & error on Mi & loss after Mi 
-) & loss after Mi 


|— )Mi +)m 2 
— )Mi +)m 2 
-)mJ+)m 2 
|— )Mi +)m 2 
h)i'i +)m 2 


| — } & error on M 2 
|+) & error on Mi 
vac) + errors on Mi and M2 
+) & errors on Mi & M 2 & loss after Mi 
|— ) & error on M 2 & loss after Mi 


|-)Mi|-)M 2 

-)mi -)m 2 

— )m 1 \— )m 2 

\-)m 1 \-}m 2 

\-}m 1 \-)m 2 



TABLE I: Possible physical scenarios corresponding to the 
measurements seen on each of the two photonic modules, Mi 
and M2. For three of the four measurement patterns, the 
correct state is listed first. 



For each case, except |— )miI — }m 2 > the correct error 
free result is listed first. Ensuring correct fault-tolerant 
operation of the perpetual network requires protecting 
the system from the following, 

1. Lost photons need to be reliably detected and re- 
placed. Additional errors may cause temporally 
correlated loss events, but these should not persist 
in the network (unless additional errors occur). 

2. Two photons should never be injected into the 
preparation network at the same time. Hence we 



tolerate a small increase in the correlated loss rate 
in order to completely suppress this possibility. 

Therefore, we always assume that no loss event has 
taken place if the modules are measured in the states 
|+)mi| — )m 2 OT |— )mi\+)m 2 an< i the photon is rerouted 
back to the start of the preparation network. The other 
possibilities, if additional errors have occurred, are, 

f . The wrong state was measured, but the photon is 
still present in the network. This is a standard 
measurement error which is effectively corrected by 
the properties of the cluster. 

2. A I vac) state was re-routed back into the network. 
This causes a temporally correlated loss event but is 
corrected (with high probability) in the next cycle. 

3. A |vac) state is re-routed back into the network 
combined with a measurement error. The topolog- 
ical cluster corrects the measurement error and a 
correlated loss error is corrected (with high proba- 
bility) in the next cycle. 

If the modules in the detection network are measured in 
the state \+) m 1 \+) m 2 or | — — )a/ 2 then we do not re- 
route back into the network and instead re-inject a new 
photon from the source (if one is available) . The adverse 
effect, when other errors occur, is that a photon may be 
unintentionally removed from the computer. Therefore, 
some measurement errors in the detection network gen- 
erate two correlated errors; the initial measurement error 
and assuming a photon is not re- injected from the source, 
the original photon is accidentally removed. 

These rules for interpreting the measurement results 
from the detection network ensures that lost photons can 
be detected and replaced, single failure events (in the 
detection modules) propagate to at most two temporally 
correlated errors and that two photons are never injected 
into the preparation network at the same time. 

C. Photon re-routing. 

Once measured by the detection network, photons are 
re-routed back into the cluster preparation network. The 
full structure of the perpetual network is illustrated in 
Fig. [8j The injection and source networks will be dis- 
cussed shortly. The return optical routing of photons 
contains a total delay of 2T. This delay is required to 
ensure that re-injected photons maintain the temporal 
arrangement of the network. Each photon spends a total 
of I0T within the cluster preparation network and the 
detection network. Once a given photon is measured by 
the second module in the detection system, it must be 
further delayed by 2T. This is due to the fact that the 
atomic system needs to be measured, taking time T, and 
to allow the classical measurement signal to be trans- 
mitted from the detector network to the source network 
(once again taking time T). The photon repetition rates 
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Cluster Preparation Network 
Total Delay = AT 

T AT 



Measurement Network 
Total Delay = 2T 




FIG. 8: Cross sectional structure of the perpetual network. The network contains three separate sub-networks. The source 
network consists of an array of photonic modules which probabilistically distill a weak coherent pulse into single photons. The 
injection network accepts photons from the source or photons recycled from the measurement network. The cluster preparation 
network deterministically couples photons and was detailed previously in Ref. [20]. Finally the measurement network which 
measures photons for computation, detects loss events and reroutes photons successfully measured back to the injection network. 



for the network are 4T and 2T respectively. Therefore, to 
maintain the temporal arrangement within the network, 
the total cycle time, Tt, must satisfy Tt mod (4T) = Tt 
mod (2T) = 0. For the network in Fig. [8] the total cycle 
time is Tt = 12T, satisfying these conditions. Conse- 
quently, photons arriving from the return path will re- 
enter the network at the correct temporal location. 



VI. SOURCE AND INJECTION NETWORK 

The final part of the perpetual design is the injection 
network. This part of the system accepts photons be- 
ing re-routed back from the measurement network and 
also accepts photons from the network of probabilistic 
sources. As explained in Sec. |IVB[ each source distills 
single photons at random, but heralded, times and the 
rest of the computer will lose photons at random, but 
heralded, times. Therefore, the injection network needs 
to be designed in a way that photons prepared by the 
sources can be routed to various injection points in the 
preparation network as they are required. 

This can be done by connecting each source to the 
optical lines immediately above. Probabilistic sources 
are therefore linked together to form a uni-directional, 
linear nearest neighbor network (uLNN). This network 
connecting sources to the preparation network is referred 
to as the shunting network. Individual wafers of pho- 
tonic modules (illustrated in Fig. p5p) have independent 



shunting networks and are not interconnected between 
separate waiters. Each of the shunting networks have a 
size along the length of each wafer related to the funda- 
mental probability that an individual source successfully 
prepares a single photon. The size of the shunting net- 
work is given by N = 0{l/p s ), where N is the number 
of sources connected in the network and p s is the proba- 
bility of distilling a single photon, per cycle. If photonic 
modules are utilized as probabilistic sources, then the 
cycle time is 3T: 2T for module initialization and mea- 
surement and an interaction time of T . The size of the 
source network is chosen such that, on average, one pho- 
ton is successfully distilled within the network every 3T. 
The source network is connected to the shunting network 
via a delay of T, this gives sufficient time for each mod- 
ule to be measured to confirm if a single photon has been 
successfully distilled before being introduced into the in- 
jection network. 

The shunting network connects each optical line to 
the one above it with an additional time delay of T. 
The boundary conditions are periodic, with the upper- 
most optical line connected to the lowermost optical line. 
Given that T is defined via the operational time of a pho- 
tonic module (between 10ns and 1/is), the spatial extent 
of the shunting network is limited to between l-100m, 
more than sufficient for a large scale array. The re- 
quired time delay of T between neighboring optical lines 
is due to the temporal arrangement of photons entering 
the preparation network, Fig. [9] illustrates. Each individ- 
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FIG. 9: Temporal arrangement of photons entering the cluster 
preparation network. Photons in neighboring optical lines are 
separated by time T. Hence if a photon is shunted one optical 
line upwards, it must also be delayed by T to be placed in the 
correct location. 

ual optical line has a pulse separation of either 2T or AT 
and each pulse in neighboring optical lines is separated 
by T. If a source produces a heralded photon at the ap- 
propriate time to be injected to its corresponding optical 
line, but the photon is not required, it is then routed 
one line up. The delay of T ensures this shunted photon 
is at the right temporal location should the photon be 
required to replace a loss event in the next optical line. 

The shunting network is designed to cycle freshly pre- 
pared photons in a loop until they are required. Each 
time a loss event is confirmed by the measurement net- 
work, the system accepts a photon from the shunting net- 
work, if one is available. For each optical line, there are 
seven distinct switching scenarios which are illustrated 
in Fig. [TUJ These switching scenarios represent all cases, 
where photons are successfully recycled, where the source 
successfully distill a state and where a photon is present 
in the shunting network. Each pattern can be briefly 
summarized, 

1. Subfigure a): Photon is successfully recycled and 
simply rerouted back into the preparation network. 

2. Subfigure b): Photon is successfully recycled and 
a second photon is present in the shunting network. 
Recycled photon is routed back into the prepara- 
tion network and the photon in the shunting net- 
work is routed up into the next optical line. 

3. Subfigure c): Photon is confirmed as lost and a 
second photon is present in the shunting network. 
The photon in the shunting network is routed into 
the preparation network. 

4. Subfigure d): Photon is confirmed as lost at the 
same time as a second photon is present in the 
shunting network and a new photon is successfully 
distilled from the source. The photon in the shunt- 
ing network is routed into the preparation network 
and the photon from the source enters the shunting 
network. 



5. Subfigure e): All three photons are present. The 
photon from the source is routed to a termination 
point and removed from the network, the other two 
photons are routed in the same way as subfigure b) . 

6. Subfigure f): Photon is confirmed as lost and 
source successfully distills a photon. This new pho- 
ton is direcly routed to the preparation network. 

7. Subfigure g): Photon is successfully recycled and 
the source distills a photon. The recycled photon is 
rerouted back into the preparation network and the 
newly distilled photon is routed into the shunting 
network. 

Photon loss can occur at any point in this network and 
we do not assume any non-demolition "probing" of the 
network to confirm if photons are still present. The only 
classical signals available are the heralding signal from 
each source and the signal from the measurement net- 
work confirming loss events. As each distilled photon is 
injected into the shunting network, they are classically 
tracked. If these photons are lost before being injected 
into the preparation network, then this will produce an- 
other loss event within the computer which, with high 
probability, will be corrected in the next cycle. 

VII. NETWORK SIMULATIONS 

To confirm that this network design operates as in- 
tended, direct numerical simulations were performed. It 
should be stressed that we are not simulating any quan- 
tum aspect of this architecture. We are simply focussing 
on the network structure under finite photon loss and low 
probability source injection. 

As detailed in previous sections, there are a total of 
12T time steps within the network, AT in the cluster 
preparation network, 2T in the measurement network, 
a AT delay between the preparation and measurement 
networks and a 2T delay when rerouting photons. During 
each of these steps, individual photons are subjected to 
loss with a probability of pl- We refer to pl as the per 
component loss probability. We also define p c as the per 
cycle loss probability, this is the probability that a photon 
is lost between entering the preparation network and re- 
entering it again after it is recycled, p c = 1 — (1 — pl) 12 ~ 

Within the source network, each photon is delayed by 
T immediately after it is distilled (in order for the module 
to be measured, confirming distillation) and again sub- 
jected to loss with probability pl- Once a photon enters 
the shunting network it is continually shunted in a loop, 
each step in the loop takes time T. The shunting con- 
tinues until the photon is accepted into the network or 
lost. 

In each simulation we vary the total time the network 
is running and the bias between p^ and the probability of 
success for each source, p s . The source probability is fixed 
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FIG. 10: Each of the seven switching scenarios dictating how photons are either recycled from the measurement network, 
accepted in or out of the shunting network or accepted from the source. In these figures the |vac) states have always been 
detected by the heralding associated with the sources or the detectors. Vacuum states can also be present which have not been 
detected. For example, a photon is lost while it is being rerouted from the detectors back to the injection network. In this case, 
the switching pattern assumes the photon is still present, re-injecting a |vac) state into the computer which will be corrected 
in the next cycle. 
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and dictated by the size of the network. For a network 
consisting of ./V optical lines, the source probability is 
p s = l/(3iV), guaranteeing that on average one photon 
is successfully distilled every 3T steps. The simulations 
are designed to determine the following, 

1. For a given network size iV = l/(3f> s ) and bias 
B = Ps/pl > 1, does the network saturate? i.e. 
how large is B such that as t — > oo each tempo- 
ral position in the network occupied with a single 
photon. 

2. For a given network size N = l/(3p s ) and a bias 
B = p s /pL, how many timesteps are required be- 
fore the network is saturated with photons. Essen- 
tially how long is required to boot-up the computer 
before computation can proceed. 

For N optical lines, the total number of photons re- 
quired to saturate the network is given by P = 9N/2. 
For a large scale computer, consisting of N x N optical 
lines, the total number of photons required to boot-up 
the computer is P — 81iV 2 /4. However, as each of the 
individual uLNN shunting networks in the computer are 
not interconnected, it is sufficient to simulate a single 
cross-section containing N optical lines. 



Network Saturation point, -9JV/2 Network + Shrrrrt Saturation, — 5iV 




Total TimeKteps, f 

FIG. 11: Simulations showing the total number of photons in 
the network as a function of total timesteps for N — 88 optical 
lines. Each curve represents higher values of 2 < B < 192 in 
increments of ten. Two bounds exists in the system. The com- 
putational network requires 9JV/2 photons to saturate the sys- 
tem. The second bound represents additional photons present 
in the shunting network. In total 5N photons are required to 
saturate the entire system. Other results for 8 < N < 168 are 
shown in Appendix. [X] 



VIII. RESULTS 

The numerical simulations are performed varying the 
network size N, the bias B and the total number of 
timesteps t. The source and loss probabilities are given 
by p s = 1/(3 AT) and p L = p s /B = 1/(3BN). Monte- 
Carlo simulations were performed using 10 3 statistical 
runs. Due to available computation power, we have sim- 
ulated up to 168 optical lines. This corresponds to a 
source probability as low as p s = 2 x 10~ 3 and hence an 
approximate error on the distilled single phot on state of 
e = 1 - F = 6.7 x 1(T 7 . Shown in Fig. [ll] is the plot 
for ./V = 88, with other simulation results shown in Ap- 
pendix. [A] This simulation illustrates the total number 
of photons present in the network as a function of t for 
2 < B < 192 averaged over 10 3 statistical runs. 

The first thing to notice is that the total number of 
photons in the network is greater than P — 9N/2 when 
B and t axe large. This is due to the shunting network. 
While the computer requires 9N/2 photons to saturate, 
additional photons are also present in the shunting net- 
work. These photons are used to replace confirmed loss 
events while the computer is in operation. On average, 
when saturated, the shunting network will contain an 
extra N/2 photons. These two bounds are illustrated on 
each plot. As B —¥ oo, the total number of photons in 
the network approaches P = 5N. 

The second set of simulations shown in Appendix. [5] 
(with Fig. 12 illustrating for A" = 88) examines the total 
number of photons present in the network as a function 
of B. Each data point is taken at the maximum value of 



t simulated in Appendix. [A] For each value of N, the net- 
work saturates once the bias reaches B 30. Note that 
this threshold bias does decrease slightly as N increases, 
but it is essentially independent of the network size. A 
bias threshold of 30 translates to a per cycle loss prob- 
ability of p c — 0Ap s . Hence, provided that the source 
probability is approximately 2.5 times higher than the 
per cycle loss probability of each photon, the system will 
eventually saturate and computation can proceed. 

Finally, we can determine an approximate boot-up 
time for the computer. For this, we fix the bias at B = 32 
and re-simulate the network for a total of 3 x 10 3 sam- 
ples. We then plot as a function of total timesteps, the 
percentage of all samples which result in complete sat- 
uration of the preparation network. Fig. |13| illustrates. 

Each curve shows a clear transition, where the total 
timesteps is large enough such that the network satu- 
rates [69 . For estimating the boot-up time for vari- 
ous network sizes, we simply find the approximate point 
where 50% of all samples lead to > 9N/2 photons in 
the network. Fig. 14 illustrates the approximate boot-up 
time for 8 < AT < 168. 

As you can see, there is a clear linear relationship be- 
tween the boot-up time and the total size of the network. 
As the size of the network is dictated by the source prob- 
ability, A^ = l/(3p s ), the time required to saturate the 
network scales inversely with p s . While, even for mod- 
estly sized networks, the total number of timesteps re- 
quired for boot-up is large, the actual physical time is 
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FIG. 12: Simulations showing the total number of photons in 
the network as a function of B for N — 88 optical lines. Each 
data point is taken at the maximum value of t simulated in 
Fig. 1 1 1 1 The system saturates at B ~ 30. Other results for 
8 < N < 168 are shown in Appendix. [B] illustrating that this 
bias threshold is essentially constant for all N. 



Network Size, N 

FIG. 14: Scaling of boot-up time as a function of network 
size, N. Each point was estimated from the simulations in 
Appendix. [A] at a bias of B = 32, finding the minimum num- 
ber of timesteps such that 50% of the simulations resulted in 
> 9N/2 photons being present in the network. The scaling 
of boot-up time scales linearly with N, hence inversely with 
p 3 . A more conservative estimate for the boot-up time for 
each datapoint results in a small constant shift in the above 
curve which, for large N, does not significantly change our 
estimates. 
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FIG. 13: Percentage of statistical samples (3x 10 3 total) which 
result in network saturation as a function of total timesteps, 
t, and network size, N. For each simulation we fix the bias at 
B — 32. Each curve shows a clear transition when the compu- 
tational network saturates and can be used for computation. 
The point at which approximately 50% of statistical samples 
result in > 97V/2 photons being present in the network is the 
boot-up time illustrated in Fig. |14| At B — 32, we are sitting 
very close to the bias threshold for low values of N. This is 
why the curves, especially for N = 8 and N = 24, show less 
than 100% of the samples saturating for large t. 



more than acceptable. Each timestep, T, represents the 
operational time for a photonic module. Hence assuming 
T = Ifis the scaling, in units of milliseconds, is approxi- 
mately, 



^boot-up (ms) 



1 



30 Ps 



(10) 



We can now summarize the requirements for a perpet- 
ual design running with highly probabilistic sources, 



Ps 


Max-p c 


^bootup (~) 


shunt network size, N, (~) 


10" 2 


4 x 10" 3 


3 ms 


30 


10" 3 


4 x 10" 4 


30ms 


3 x 10 2 


10~ 4 


4 x 10~ 5 


300ms 


3 x 10 3 


10~ 5 


4 x 10" 6 


3s 


3 x 10 4 


10" 6 


4 x 10" 7 


30s 


3 x 10 5 



TABLE II: Estimation of boot-up times and shunting network 
sizes for various values of distillation probability, p a . Max-p c 
is the maximum per cycle probability of photon loss such that 
the computational network will saturate (corresponding to a 
Bias of approximately p s = 2.5p c . 



A reasonable parameter range for topological cluster 
states requires a per cycle loss probability, p c , between 
approximately 1CP 2 and 10 -4 . Higher loss probabilities 
begin to require unreasonable cluster resources to cor- 
rect [S7j while loss probabilities lower than m 10 -4 are 
experimentally unrealistic. Therefore, source probabili- 
ties can be in the range of approximately 2.5% — > 0.025% 
with each shunting network connecting 10 — > 300 optical 
lines and an boot-up time of between 1 — > 100ms. 

These estimates represent approximate upper bounds 
for the boot-up time and the size of each shunting net- 
work. If each probabilistic source succeeds with a much 
higher probability (while loss probabilities in the network 
remain fixed), then both the boot-up time and the size 
of each shunting network will decrease. 
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IX. CONCLUSIONS 

This discussion has shown how highly probabilistic 
photon sources can be integrated into a scalable opti- 
cal architecture. By replacing deterministic sources and 
detectors with the photonic module, this modified design 
is constructed exclusively from a single quantum compo- 
nent. These results equally apply to any other form of 
highly probabilistic photonic sources, provided they can 
be integrated into a computational network similar to the 
optical architecture used in this work. 

The introduction of the shunting network, connecting 
probabilistic sources to the cluster preparation network, 
combined with the ability to recycle measured photons 
essentially leads to a pseudo set of on demand sources. 
It is easily shown that by having sources with a slightly 
higher probability of success than loss, the system will 
be able to compensate photon loss during computation. 
The shunting network solved the problem of being able 
to route photons prepared in random locations in the 
network to other random locations where photons are 
lost. This structure effectively achieved the same goal as 



extremely complicated multi-port switches and extensive 
optical delay. 

This analysis demonstrates the practicality of a lim- 
ited amount of highly probabilistic quantum components 
in a large scale architecture. Unlike other results, illus- 
trating the theoretical ability to construct a quantum 
computer from probabilistic components, this work has 
maintained the explicit architectural structure of the op- 
tical computer. By starting from a fully deterministic ar- 
chitecture, we could carefully integrate a limited amount 
of probabilistic technology without sacrificing the overall 
structure of the system. The next major step is to de- 
termine if additional probabilistic technology can replace 
other currently deterministic components. The results 
presented in this paper give us a certain amount of opti- 
mism that this can be achieved. 
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As we set B = 32, we are sitting very close to the bias 
threshold for small values of N , hence for large t not 
all samples result in network saturation. This is why the 
curves, especially for N = 8 and N = 24, show less than 
100% saturation for large t 
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APPENDIX A: FURTHER SIMULATIONS: 
NUMBER OF PHOTONS AS A FUNCTION OF 
TIMESTEPS 



Shown here are further simulations for 8 < N < 168, 
showing the total number of photons in the network as 
a function of total timesteps, t, and bias, B. Total sam- 
ples are 10 3 and for larger networks the maximum bias 
simulated is reduced due to computational resources. 



APPENDIX B: FURTHER SIMULATIONS: 
NUMBER OF PHOTONS AS A FUNCTION OF 
BIAS 



Shown here are further simulations for 8 < N < 168, 
showing the total number of photons in the network as a 
function of bias, B. Each data point is taken at the max- 
imum value of t simulated in Appendix. |Aj For each net- 
work size, the system saturates at B ss 30. The threshold 
bias decreases slightly as N increases. 
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FIG. 15: N = 8. Simulations of 2 < B < 192. As the total 
number of photons is quite low, this simulation suffers from 
more statistical variance for 10 3 samples. 
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FIG. 16: N = 24. Simulations of 2 < B < 192 
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FIG. 17: N = 40. Simulations of 2 < B < 192. 
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FIG. 18: A^ = 56. Due to computational resources, we restrict 
simulations to 2 < B < 172. 



FIG. 21: A" = 120. Due to computational resources, we re- 
strict simulations to 2 < B < 132. 
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FIG. 19: iV = 72. Due to computational resources, we restrict 
simulations to 2 < B < 132. 
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FIG. 22: A^ = 136. Due to computational resources, we re- 
strict simulations to 2 < B < 102. 
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FIG. 20: A^ = 104. Due to computational resources, we re- 
strict simulations to 2 < B < 132. 



FIG. 23: A" = 152. Due to computational resources, we re- 
strict simulations to 2 < B < 92. 
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FIG. 24: N = 168. Due to computational resources, we re- 
strict simulations to 2 < B < 82. 
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FIG. 27: AT = 40. Simulations of 2 < B < 192. 
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FIG. 25: N = 8. Simulations of 2 < S < 192. 
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FIG. 28: AT = 56. Due to computational resources, we restrict 
simulations to 2 < B < 172. 
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FIG. 26: JV = 24. Simulations of 2 < B < 192. 
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FIG. 29: iV = 72. Due to computational resources, we restrict 
simulations to 2 < B < 132. 
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FIG. 30: N = 104. Due to computational resources, we re- 
strict simulations to 2 < B < 132. 
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FIG. 33: N = 152. Due to computational resources, we re- 
strict simulations to 2 < B < 92. 
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FIG. 31: iV = 120. Due to computational resources, we re- 
strict simulations to 2 < B < 132. 
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FIG. 34: N = 168. Due to computational resources, we re- 
strict simulations to 2 < B < 82. 
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FIG. 32: N — 136. Due to computational resources, we re- 
strict simulations to 2 < B < 102. 



